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From both theoretical and experimental points of view symmetric states constitute an important class of multi- 
partite states. Still, entanglement properties of these states, in particular those with positive partial transposition 
(PPT), lack a systematic study. Aiming at filling in this gap, we have recently affirmatively answered the open 
question of existence of four-qubit entangled symmetric states with positive partial transposition and thoroughly 
characterized entanglement properties of such states [J. Tura et al, Phys. Rev. A 85, 060302(R) (2012)] With 
the present contribution we continue on characterizing PPT entangled symmetric states. On the one hand, we 
present all the results of our previous work in a detailed way. On the other hand, we generalize them to sys- 
tems consisting of arbitrary number of qubits. In particular, we provide criteria for separability of such states 
formulated in terms of their ranks. Interestingly, for most of the cases, the symmetric states are either separable 
or typically separable. Then, edge states in these systems are studied, showing in particular that to character- 
ize generic PPT entangled states with four and five qubits, it is enough to study only those that assume few 
(respectively, two and three) specific configurations of ranks. Finally, we numerically search for extremal PPT 
entangled states in such systems consisting of up to 23 qubits. One can clearly notice regularity behind the 
ranks of such extremal states, and, in particular, for systems composed of odd number of qubits we find a single 
configuration of ranks for which there are extremal states. 



I. INTRODUCTION 

Characterization of entanglement [ 1 ] in composite quantum 
states with positive partial transposition (PPT states) remains 
a difficult problem. One of the reasons for that is the lack of 
a universal separability criterion allowing to distinguish un- 
ambiguously separable from PPT entangled states (see, nev- 
ertheless, e.g. Ref. [2] for numerous necessary separability 
conditions). There are, however, methods providing some in- 
sight into the structure of PPT entangled states. One of them 
exploits the fact that all states that remain positive under par- 
tial transposition form a convex set, which as a proper subset 
contains the PPT entangled states. To fully characterize the 
latter, it is then enough to know all the extremal points of this 
convex set. This approach has recently been extensively stud- 
ied (see Refs. [3-7]). In particular, it allowed to solve the open 
problem of existence of four-qubit PPT entangled symmetric 
states [8], and also, although in an indirect way, disprove the 
Peres conjecture in the multipartite case [9]. 

The problem of characterization of PPT entangled states 
is even more complicated in the multipartite case. Clearly, 
the set of PPT states arises by intersecting sets of states that 
remain positive under partial transpositions with respect to 
single bipartitions, thus its boundary becomes more compli- 
cated with the increasing number of parties. Nevertheless, 
the complexity can be reduced by imposing some symme- 
tries. For instance, demanding that the states under study com- 
mute with multilateral action of unitary or orthogonal groups 
leads to classes of multipartite states whose full characteriza- 
tion with respect to entanglement becomes possible (see e.g. 
Refs. [10]). 

Another interesting example of a class of states obtained 
by imposing some symmetry are those supported on the sym- 
metric subspace of a given multipartite Hilbert space. The 
so-called symmetric states have recently been attracting much 
attention [11-18]. In particular, the underlying symmetry 



allowed for the use of the Majorana representation [19] for 
an identification of SLOCC classes of multipartite symmetric 
states [14] (see also Refs. [15]). The same symmetry provides 
advantages in calculating certain entanglement measures [16]. 
Another motivation comes from the recent experimental real- 
izations of symmetric states of many-qubits, as for instance, 
the six-qubit Dicke states [20] or the eight-qubit GHZ states 
[21] (see also Ref. [22] in this context). 

However, more effort has been devoted to the pure sym- 
metric states, leaving the characterization of entanglement of 
mixed, in particular PPT states as an open problem. It is 
known so far that for N = 2, 3 all PPT symmetric states are 
separable [11]. Then, examples of five or six-qubit PPT entan- 
gled symmetric states were found in Refs. [17, 18]. Recently, 
the remaining case of A" = 4 has been studied in Ref. [8], 
where the open question as to whether partial transposition 
serves in this case as a necessary and sufficient condition for 
separability (as this is the case for A^ = 2, 3) has been given a 
negative answer. The main aim of the present paper is to con- 
tinue the characterization of PPT entanglement in symmetric 
states. We discuss in detail methods used in Ref. [8] and then 
generalize them to the case of arbitrary N. We derive separa- 
bility criteria for PPT symmetric states in terms of their ranks 
and ranks of their partial transpositions. Then we exclude con- 
figurations of ranks for which they are genetically not edge. 
Finally, we adapt to the multipartite case an algorithm allow- 
ing to search for extremal PPT entangled states [3] (see also 
Ref. [7]). Exploiting it, we study ranks of the extremal PPT 
entangled symmetric states consisting of up to 23 qubits. In- 
terestingly, we show that there are at most three distinct con- 
figurations of ranks for which we find extremal PPT entangled 
symmetric states, and, in particular, for odd A~ there is only a 
single such configuration. 

The paper is structured as follows. In the next section (Sec. 
II) we recall all notions and facts necessary for further consid- 
erations. In Sec. Ill we investigate the entanglement proper- 
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ties of the PPT symmetric states. Then, in Sec. IV we seek 
extremal entangled PPT symmetric states consisting of even 
more than 20 qubits and classify them with respect to their 
ranks. The results obtained for exemplary systems consisting 
of four, five and six qubits are collected in Sec. V. We con- 
clude in Sec. VI. 



II. PRELIMINARIES 

Let us start with a couple of definitions that we will use 
throughout the paper. Let 

H N = € dl ®...®<D dN (1) 

denote a multipartite product Hilbert space and D the con- 
vex set of density operators acting on Hn- By R{p), r (p), 
K(p), and k(p) we denote, respectively, the range, rank, ker- 
nel, and the dimension of the kernel of a given p G D. Then, 
A\ , , . , , An will stand for the subsystems of a given TV-partite 
p, and, in the case of low N, we will also denote them by 
A, B, etc. 

PPT and separable states. Let us now split the set / = 
{Ai, . . . , An} into two disjoint subsets S and S (S U S = I) 
and call it bipartition S|S. We say that a given state p acting 
on Hn is PPT with respect to this bipartition iff p Ts > 0. 
Clearly, states with this property make a convex set denoted 
Ds- An element of D whose partial transpositions with re- 
spect to all bipartitions (notice that for a given bipartition S|S, 
partial transpositions with respect to S and S are equivalent 
under the full transposition) are positive will be called fully 
PPT, and, since in this paper we deal only with fully PPT 
states, we will be calling them simply PPT states. Clearly, 
such states make also a convex set which is simply the inter- 
section of Ds for all S. 

A particular example of a state that is PPT is the fully sepa- 
rable state [23, 24]: 

P = Yl Pi f ) A 1 ®-..® Pa n , Pi>0, $3pi = l, (2) 

i i 

where p\, denote density matrices representing all subsys- 
tems. Clearly, in multipartite systems one may define vari- 
ous types of separability (see e.g. [24, 25]). Nonetheless, as 
we will see later, in the symmetric case a given state p is ei- 
ther genuine multipartite entangled, i.e., cannot be written as a 
convex combination of states which are separable with respect 
to in general different bipartitions, or takes the form (2). 

Edge states. An important class of entangled PPT states are 
the so-called edge states [26-28]. We call p acting on Hn 
edge iff there does not exist a product vector |ei) <8> . . . 8> |ejv) 
with je») G € d ' such that |ei)(g>. . .<g>\e N ) G R(p) and (|ei)(g> 
|eAr)) Cs G R{p Ts ) for all S, where by Cs we denoted 
partial conjugation with respect to S. The importance of edge 
states in the separability problem comes from the fact that any 
PPT state can be decomposed as a mixture of a fully separable 
and an edge state [26]. Alternatively speaking these are states 
from which no product vector can be subtracted without losing 
the PPT or positivity property, meaning that they lay on the 



boundary of the set of PPT states. However, they do not have 
to be extremal, although any extremal state is also edge. 

Edge states have been studied in bipartite or three-partite 
systems and many examples have been found [27-29]. 

Symmetric states. Let us now concentrate on the iV-qubit 
Hilbert space 

n 2 , N = {<d 2 )® n (3) 

and consider its subspace Sn spanned by the unnormalized 
vectors 

\E?) = |{0, i},{l, N-i}) (i = 0,...,N), (4) 

which are just symmetric sums of vectors being products of 
i zeros and N — i ones. These vectors, when normalized, 
are also known as Dicke states. For further benefits, let us 
notice that the dimension of Sn is N + 1, and therefore it is 
isomorphic to <C N+1 , which we denote Sn — (D +1 . Also, 
by Vn we will be denoting the projector onto Sn- 

We call a state p acting on H^n symmetric iff it is sup- 
ported on Sn, or, in other words, R(p) C Sn- In yet another 
words, p is symmetric iff the equations 

V aP = pVl = P (5) 

are obeyed for any permutations a, a' G Sat, where Ejv is the 
group of all permutations of an iV-element set, while V a is an 
operator defined as V^i) ■ • ■ \iPn) = |VV(i)} ■ ■ • IVV(iV)) for 
any vectors G C 2 . 

In the case of symmetric states the number of relevant par- 
tial transpositions defining the set of PPT symmetric states 
£>p V pT is significantly reduced. This is because positivity of 
a partial transposition with respect to some subset S is equiv- 
alent to positivity of all partial transpositions with respect to 
subsystems of the same size \S\. Together with the fact that 
for a given bipartition S\S, p Ts > <^=> p T ~ > 0, one has 
[TV/2J partial transpositions defining Dpp^,. We choose them 
to be Ta 1 , Ta 1 a 2 , e t c -> however, for simplicity we will also be 
denoting them as T\ = Ta x , T-2 = Ta 1 a 2 , an d so on. Alterna- 
tively, in systems of small size, we will use Ta, Tab, Tabc, 
etc., to denote the relevant partial transpositions. 

Let us now notice that since Si = (D I+1 , a iV-qubit sym- 
metric state p can be seen with respect to a bipartition S\S as 
a bipartite state acting on (D' s ' +1 (g> (D JV_ I S I +1 . This gives us 
nontrivial bounds on the ranks of partial transpositions with 
respect to all S, namely, 

r(p Ts )<(\S\ + l)(N-\S\ + l) (6) 

for | SI = 0, . . . , [N/2\, which, in particular means that 
r (p) < N + 1. A very convenient way of classifying PPT 
states is through their ranks and ranks of their partial transpo- 
sitions, i.e., the |_JV/2j -tuples 

(r(p), r(p T -^ ),r(p TA ^), r^ 1 - 1 ^ )) 

= (r(p), r(p T i),r(p T * ),..., r(p T ^)) . (7) 

Finally, let us recall that p acting on some bipartite Hilbert 
space % = < C dl <8> <D rf2 is said to be supported on iff 
R(Pa) = C^ 1 and R{ps) = <C d2 - Alternatively speaking, p 
is not supported on TL2 if either pa or ps has a vector in the 
kernel. 
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III. CHARACTERIZING PPT ENTANGLEMENT IN 
iV-QUBIT SYMMETRIC STATES 

Here, exploiting the results of Refs. [30, 31] we derive sep- 
arability criteria for the PPT symmetric states in terms of the 
ranks (7). Recall that in these papers it was shown that any 
PPT state p supported on a Hilbert space <C dl <Z> <C d2 is sepa- 
rable if r(p) < max{<ii, e^}- 

Then, we study the edge symmetric states and in particular 
we show that symmetric states assuming certain ranks cannot 
be edge. 

A. Separability 

To begin with the separability properties let us recall that 
if a pure iV-partite symmetric state \tp) is separable with re- 
spect to some bipartition, then it must be fully separable (2), 
i.e., IV-} = \e)® N with |e) G € 2 (see Refs. [11, 13]). This 
straightforwardly implies that entangled symmetric pure and 
thus mixed states have genuine multipartite entanglement (see 
also Ref. [13]). Indeed, if a symmetric p can be written 
as a convex combination of density matrices, each separable 
across some, in general different, bipartition, then p has pure 
separable vectors in its range. Since each such vector is sym- 
metric, it assumes the above form \e)® N , meaning that p is 
fully separable. Thus, throughout the paper, by saying that a 
symmetric state p is separable we mean that it is fully separa- 
ble, i.e., it takes the form (2). 

Let us now establish some conditions for separability in 
terms of ranks of p. We start with the following technical 
lemma. 

Lemma 1. Consider an N-qubit symmetric state and a bipar- 
tition S\S (\S\ < N — \S\). Then, let ks and kg (r$ and r-g) 
denote the dimensions of the kernels ( ranges) of subsystems of 
p with respect to S\S. The following statements hold: 

• ifk-g>0 then r(p) < r-g, 

• ifks>0 then r(p) < r,s, 

• ifks > and kg > 0, then r(p) < min{rs, r-g}. 

Proof. Recall first that with respect to the bipartition S\S, the 
iV-qubit symmetric state p can be seen as a bipartite state act- 
ing on <Dl s l +1 ® = Cl s l +1 ® <C N ~\s\+\ 

We will prove the first case and then the remaining two will 
follow. Assume then that kg- > 0, meaning that p-g has k-g 
linearly independent vectors \cf>i) (i — 1, . . . , kg) in the ker- 
nel. Consequently, for any |5'|-qubit symmetric vector the 
projected vectors V N (\i>)\<f>i)) (i = 1, . . .,N-r§- \S\ + 1) 
belong to K(p). In what follows we will show that by choos- 
ing properly vectors \ip), one is able to find N — r-g + 1 linearly 
independent vectors in K(p) of this form. 

First, we prove that for any , N — r-g — \S\ + 1 pro- 
jected vectors VN{\ip)\<t>i)) G K{p) are linearly indepen- 
dent. Towards this end, let us assume, in contrary, that there 
exists a collection of nonzero numbers on G (D such that 



J2i a iPN{\i>) \4>i)) = 0- The latter is equivalent to saying 
that the vector \ip) ® ^ a.i\4>i) sits m the kernel of Vn- But 
since oti\4>i) is an IS'l-qubit symmetric vector, this is pos- 
sible only if are linearly dependent, contradicting the fact 
that they span the kernel of Ti'sP- 
Now, we consider particular vectors 

W)=V N (\Ef l )\ct> i )) (8) 

with j = 0, . . . , |5| and i = 1, . . . , \~3\ - + 1. As al- 
ready proven, for any j, the vectors |$^) make an (IS"! — 
r-g + l)-element linearly independent set. Let us now con- 
centrate on the vectors |<£>°) and choose the one for which 
(E^||<I>°) is nonzero, say |$^). Notice that by the construc- 
tion (E^\$°) = for k > \S\; to obtain |$?) we symmetrize 

\<j)i) with \E^) = |0}®l s l, meaning that |$°) decomposes 
into the symmetric vectors l-E^} with k < \S\. If, how- 
ever, (_E^|<£>°) = for all i, we choose the one for which 

{E^^_ 1 ^ 0, etc. Clearly, repeating this we must find the 

desired vector as otherwise |<I>"} = for all i, contradicting 
the fact that the vectors |<£>°) are linearly independent. 

Let us then assume for simplicity that |$°) is such that 
(E*L \$<?) ^ 0, meaning that ^ [cf. Eq. (8)]. 

Consequently, the vectors (j = 1, . . . , 15*1), when de- 
composed into the symmetric basis of Sn, contain \E^) with 
k > \S\ + 1 and therefore are linearly independent of the set 
Moreover, by the very construction, they make an 
\S | -element set of linearly independent vectors themselves, 
meaning that the vectors |$°) (i = 1, . . . , \S\ — r-g + 1) 
together with | <I>^ } (j = 1 , . . . , | S | ) make the desired set of 
TV — r-g + 1 linearly independent vectors in K(p). 

Consequently, k(p) > N — r-g + 1 which, taking into 
account the maximal possible rank of a symmetric state p, 
gives the bound r(p) < N + 1 — N + r-g — 1 = r-g. In 
an analogous way one proves the second case, i.e., when 
ks > 0. Precisely, following the above arguments, one sees 
that kg linearly independent vectors in the kernel of pg gives 
at least TV — r$ + 1 linearly independent vectors in K(p), 
imposing the bound r(p) < r$. To prove the third case, 
one just chooses the tighter of both the above bounds, i.e., 
f(.P) < rninjrs, r-g}. □ 

Essentially, this lemma says that if the symmetric state p 
is not supported on (D' s ' +1 ® (C^l 5 ^ 1 with respect to the 
bipartition S\S, its rank is bounded from above by ranks of 
its subsystems. In the particular case when the subsystem S 
consists of a single party, it straightforwardly implies that if 
r (p) > N then p has to be supported on <C 2 ®<C N with respect 
to the bipartition one versus the rest (Ai \A2 . . . An). 

The following fact was already stated in Ref. [11], however, 
a detailed proof was not given. We will exploit lemma 1 to 
demonstrate it rigorously. 

Theorem 1. Let p be a N-qubit PPT symmetric state. If it is 
entangled then r(p) = N + 1, i.e., p is of maximal rank. 
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Proof. An iV-qubit symmetric state can be seen as a bipartite 
state acting on (D 2 ® <C N with respect to the bipartition one 
qubit versus the rest, as for instance ^4i|^2 . . . An. Let us 
denote by pa x and pa x ...a n the subsystems of p with respect 
to this bipartition. Assuming then that r{p) < N, the results 
of Ref. [30] imply that p is separable provided it is supported 
on (D 2 (g> <D . If, however, the latter does not hold, there are 
vectors in the kernel of either pa 1 or pa 2 ...a n - In the first 
case, lemma 1 implies that r(p) = 1 and p is a pure product 
vector, while in the second case r(p) is upped bounded by 
the rank of pa 2 ...a n - Again, results of Ref. [30] apply here, 
meaning that p is separable. □ 

Alternatively speaking, this theorem means that there are 
no PPT entangled A-qubit symmetric states of rank less than 
N + 1. On the other hand, it provides only a sufficient con- 
dition for separability, as there are separable symmetric states 
of rank N + 1. Another consequence of theorem 1 is that for 
an arbitrary bipartition S\S, a PPT entangled state p and its 
partial transposition p Ts are supported on the corresponding 
Hilbert space (D |s|+1 (g)(D w ~ |5|+1 . Specifically, one has the 
following lemma. 

Lemma 2. Consider an N -qubit PPT symmetric state and 
an arbitrary bipartition S\S. If p is entangled, then p Ts is 
supported on the bipartite Hilbert space corresponding to the 
bipartition S\S, i.e., (TJ^ I+1 ® <C N ~^ +1 . In other words, if 
for some bipartition S\S, p Ts is not supported on (D' s ' +1 ® 
<C N ' S ' +1 , then p is separable. 

Proof. Assume that the PPT state p is entangled but the partial 
transposition p Ts is not supported on the Hilbert space corre- 
sponding to the bipartition S\S, i.e., (D^+^C^- 1 ^^. This 
means that one of its subsystems, say the S one, contains an 
S* | -qubit symmetric vector \(j>) in the kernel. Consequently, 
for any |5|-qubit vector \ip), the following implication 

p Ts \i>)\<t>)=o P \rm=o (9) 

holds. Putting for instance \ip) = |0}®' s ', one sees that the 
symmetrized vector 'Pjv(|0)®'' s ' <g> \4>)) belongs to K{p). As 
a result r(p) < N and theorem 1 implies that p is separable, 
leading to the contradiction. □ 

Then, with the aid of lemma 2, we can prove the analog of 
theorem 1 for the ranks of partial transpositions of p. 

Theorem 2. Let us consider an N -qubit PPT symmetric state 
p and a bipartition S\S with an arbitrary S (\S\ < \S\). If it 
is entangled then r{p Ts ) > N — \S\ + 1. In particular, if p is 
entangled then r(p TA ) > N + 1. 

Proof. Due to lemma 2, we can assume that with respect to 
the bipartition on S\S, p is supported on (D |5|+1 (g> (D Ar ^l s l+ 1 
as otherwise it is separable. Then, if r(p Ts ) < N — \S\ + 1, 
the results of Ref. [31] imply that p is separable. Noting that 
this reasoning is independent of the bipartition, we complete 
the proof. □ 

In other words, any PPT symmetric states whose rank obeys 
r(p Ts ) < N — \S\ + 1 for some bipartition, is separable. 



Still, by using more tricky bipartitions we can provide fur- 
ther separability conditions for generic symmetric states in 
terms of the ranks. In this direction we prove the following 
theorem. 

Theorem 3. Consider an N -qubit PPT symmetric state p and 
a bipartition S\S. Ifr{p Ts ) < (\S\ + 1)(N - \S\), then the 
generic p is separable. 

Proof. Consider an A-qubit state a = p Ts which is PPT, but 
no longer symmetric, and a party which does not belong to S, 
say An- With respect to the bipartition An versus the rest, 
a can be seen as a bipartite state acting on (Tj(l' s l +1 )' Ar- l' s l) <g) 
(D 2 . Since p is fully PPT, it clearly follows that ct Ta « > 0. 
This, together with the fact that r(p Ts ) < (\S\ + 1)(N - 
\S\), implies that a has to be separable across the bipartition 
A x . . . An-^An, i.e., 

& = p Ts = IV'iX^iUx—AAr-i ® \ei){ei\ AN , (10) 

i 

provided that it is supported on (D(l s l+ 1 )( Ar ~l s 'l) <g> (D 2 , which 
generically is the case. 

Now, let us notice that p Ts is still symmetric with re- 
spect to subsystem S, i.e., V-gp Ts V-g = p Ts . Since then 
all vectors |V , j)l e i) appearing in the decomposition (10) be- 
long to R(p Ts ), they also enjoy the above symmetry, i.e., 
r Ps\ei)A N \i>i)A 1 ...A N - 1 = \e i )A N \4>i)A 1 ...A N - 1 for all i. 
This, as it is shown below, implies that for any i, \ipi) = 

|V^)l e i)® with IV'i) being from the Hilbert space corre- 
sponding to the subsystem S. Putting these forms to Eq. (10), 
one arrives at 

p Ts =J2^ i \ s ^>\ e ^ e ^ 1 - 

i 

Now, one can move the partial transposition with respect to S 
to the right-hand side of the above identity, and use the fact 
that p is symmetric, which leads us to the form (2). 

To complete the proof let us show that if 

~Pa\jsW)a\^)siss = \ e )A\^)sus holds for some one " 
qubit and A^-qubit vectors |e) and \ip), then the latter assumes 

the form \4>}g uS = |e) (8, l' s 'l with \ip) belonging to the 

Hilbert space associated to the subsystem S. Here by Vaus 

we denote a projector onto the symmetric subspace of the 

Hilbert space corresponding to the qubits A and S. 

To this end, let us decompose 

W) = Y / \ i h\ < t , i)s, d2) 

i 

where {\i)} denotes any orthogonal basis in the Hilbert 
space associated to the subsystem S. It is clear 
that VausI^aI^sus = \e)A\tp)sus is equivalent to 
VAus\e)A\<Pi)s = \e)A\<f>i)s for all i. By virtue of the results 
of Refs. [11, 13], the latter can hold only if = |e)®l s l 
for any i. For completeness, let us recall the proof of this fact. 
For this purpose, assume that |e)|</>) e £|s|+i for some one- 
qubit and 1 5 1 -qubit vectors \e) and \<j>), respectively. Then, 



5 



one immediately concludes that \<j>) must be symmetric and 
consequently can be written as 



ISI+l 



= E 



\S\\ 



(13) 



with a,i G (D. Putting |e) = (a, 6) and utilizing the fact that 

P| S | + i|e)|0) = |e> |0>, one has 



|S|+i 
E a * 



i=l L V i-l 



(J-l) iplSI+lv , (» |5 l) | F |S|+1\ 
ISI+1N I- * ' + (\S\+1\ 1 ' 



•Sl + l 



(O|0)+6|l))® J2 a *\ E : 



\s\- 



(14) 



Projection of the above onto vectors |0)®d s 'l +1 ), |0. . .01), 
|0 . . . Oil), etc., leads to equations aacj = bctj-i (j = 
2, . . . ,\S\ + 1), which, in turn, imply that oij = (6/o) J ' — 1 ai 
C? = 2,...,|5| + l). 
Putting now | </),;) = |e)®l s for all i to Eq. (12), one gets 

U') = E i^i^}s = E = w<rie>® |s| , (is) 

which finishes the proof. □ 

Notice that the above theorems imply that for most of 
the possible ranks, the symmetric states are either separa- 
ble or genetically separable. The maximal rank of a par- 
tial transposition with respect to a given subsystem S is 
(\S\ + 1)(N - \S\ + 1). Taking into account the fact that 
we can put r(p) = N + 1 (as otherwise the state is separable, 
cf. theorem 1), we have in total 



[N/2\ 



[ (|S| + l)(/v-|S| 



7V!(iV/2 + l) 2(LAr/2J - Ar/2)+1 



|S| = 1 

(16) 

possible configurations of the relevant ranks (7) that can be 
assumed by the symmetric states. With theorems 2 and 3, 
we see that for a transposition with respect to S, symmet- 
ric states of the first (\S\ + 1)(N — \S\) of the correspond- 
ing ranks are either separable or generically separable. This 
leaves only \S\ + 1 of ranks with respect to this bipartition for 
which they do not have to be generically separable. Taking 
into account all the relevant partial transpositions, we have 
in total ([N/2\ + 1)1 of the remaining cases where one can 
search for PPT entangled symmetric states. This is clearly a 
small portion (rapidly vanishing for large N) of all the possi- 
ble ranks [cf. Eq. (16)]. For instance, for N = 4,5, this gives 
us treatable 6 configurations of ranks (out of all, respectively, 
72 and 120 obtained from Eq. (16)), while for TV = 6, 7 this 
number amounts to 24 (Eq. (16) gives in these cases, respec- 
tively, 2880 and 5040 possible ranks). As we will see later, if 
additionally we ask about possibility of being an edge state, 
these numbers may be further reduced. 



B. TV-qubit symmetric edge states 

Let us now single out the configurations of ranks where the 
symmetric states can be edge. Clearly, theorem 3 implies that 
the ranks of p Ts for aU S must be larger than (|S'| + 1)(/V-|S'|) 
as otherwise the generic symmetric states are separable. In 
what follows we will provide a few results allowing to bound 
the ranks from above. 

In general, as already discussed in Sec. II, to prove that a 
given p is not edge one has to prove that there is |e) G (D 2 
such that \e)® N G R(p) and [\e)® N ] Cs G R{p Ts ) for all S. 
Assuming that the rank of p is maximal, r(p) = N + 1, the 
above is equivalent to solving of a system of M/° Ts ) 
equations 



[(e|n Cs |*f)=0, 



(17) 



where |*f) G K(p Ts ) and S = A,AB,.... By putting 
|e) = (1, a), one reduces Eqs. (17) to a system of polynomial 
equations P(a, a*) — in a and a* . This is clearly a hard 
problem to solve (see, e.g., the discussion in Ref. [30]). Still, 
under some assumptions and using a method of Ref. [33], it 
is possible to find a solution to a single equation of that type. 



Lemma 3. Consider an equation 



k 

E 

i=0 



(a*yQ l (a)=0 



(18) 



where Qi (i = 0, . . . , k) are some polynomials. If 
maxijdegQi} = degQ^, = n > k and degQo ~ m > k, 
then this equation has at least one solution. 

Proof. Notice that, via the results of Ref. [30], Eq. (18) has 
generically at most 2 fc ~ 1 [k + n(n — k + 1)] complex solutions. 
To find one, in Eq. (18) we substitute a = rs and a* = r/s 
with r G R and s G <D, obtaining 



^/"V«Q,(r S )=0. 



(19) 



i=0 



Treating r as a parameter and s as a variable, our aim now is to 
prove that for some r there is s such that \ s\ = 1 and Eqs. (18) 
and (19) is obeyed. For this purpose, let us first put s = x/r 
with x G <D, which gives us 



E*" 



V 2i Qi(i) = 0. 



(20) 



In the limit r — > oo, the left-hand side of the above equa- 
tion approaches Qk{x), meaning that Eq. (19) has n solutions 
s i° ~ ► (* = 1, . • . , n). Then, in the limit of r — > 0, the 
left-hand side of Eq. (20) goes to Qo(a), implying that Eq. 
(19) has m solutions — > oo (i = 1, . . . , m). 

Then, one sees that Eq. (19) has at most n + k solutions 
with respect to s. Consequently, for r — > oo and r — » 0, 
Eq. (19) has additional k roots s°° (i = n+ 1, . . . , n + k) and 



-k — m roots s? (i = m+ 1, . 



-k), respectively, which 
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can remain unspecified. As r varies continuously from zero to 
large values, all roots s° must continuously tend to s°°. But, 
since m > k, at least one of m roots — >• oo (i = 1, . . . , m) 
must tend to one of the n roots s°° (i = 1 , . . . , n) which are 
close to zero. This means that there is at least one pair (r, s) 
with \s\ = 1 solving Eq. (19) and thus Eq, (18). □ 

With the aid of the above lemma we can prove the following 
theorem. 

Theorem 4. Let us consider N-qubit PPT symmetric state 
and a subsystem S of size 1 < \S\ < \N/2~\ — 1, and assume 
thatr(p Tx ) = (\X\ + 1)(N- \X\ + 1) (maximal) for all X 
except for X = S for which r(p Ts ) = (\S\ + l)(N - \S\ + 
1) — 1. Then, generically such states are not edge. 

Proof. We prove that under the above assumptions it is gener- 
ically possible to find a product vector |e)® Ar G R(p) such 
that (\e)® N ) Cx G R{p Tx ) for all subsystems X. Clearly, 
all ranks of p are maximal except for the one corresponding 
to the partial transposition with respect to the subsystem S, 
which is r{p Ts ) = (\S\ + l)(iV - \S\ + 1) - 1 (maximal di- 
minished by one). Denoting by |^>) the unique vector from 
the kernel of p Ts , one then has to solve a single equation 
(#|(| e *)®|S| $ | e )®(W-|s|)) = o After putting |e) = (l,a) 
with tt£(C, the latter can be rewritten as 
|S| 

J2(a*TQi(a) = 0, (21) 

i=0 

where Qi(a) (i = 0, . . . , \S\) are polynomials of degree at 
most N — \S\ and generically they are exactly of degree N — 
\S\. Due to the assumption that N — \S\ > \S\, lemma 3 
applies here, implying that (21) has at least one solution and 
generic symmetric p of the above ranks is not edge. □ 

The above theorem says that generic PPT symmetric states 
having all ranks maximal except for a single one correspond- 
ing to a partial transposition with respect to a subsystem S 
such that 15*1 < N — \S\, for which the rank is one but max- 
imal, are not edge. This method, however, does not work for 
even N and in the case when the chosen partial transposition 
is taken with respect to the half of the whole system, i.e., 
I^l = 15*1 = N — | S\. This is because in this case the re- 
sulting equation (20) is of the same orders in a and a* and 
the above method does not apply. 

Still, however, using a different approach, we can prove an 
analogous fact for N = 4 and N = 6. Specifically, we will 
show that all symmetric states of all ranks maximal except for 
the last one, which is one less than maximal, are not edge. To- 
wards this end, let us start with some general considerations. 

Let p be an iV-qubit symmetric state with even N, and let 
now S denote a particular bipartition consisting of first N/2 
qubits (half of the state). Assume then that all ranks of p are 
maximal except for the one corresponding to the partial trans- 
position with respect to S, for which it is one but maximal, 
i.e., (N/2 + l) 2 — 1. Consequently, there is a single vector 
\^} in the kernel of p Ts , meaning that p is not edge iff a single 
equation [cf. Eqs. (17)]: 



with |e) G (D 2 has a solution. For this purpose, one notices 
that 

V s,sP TsV s,s = ^ Ts > (23) 

with V s g denoting a unitary operator swapping the subsys- 
tems S and S (notice that both are of the same size). Then, 
because is the unique vector in K(p Ts ), it has to enjoy the 
same symmetry, i.e., V s 5!^) = |^*)- Denoting now by My 
the matrix representing elements of in the product basis 
\E^ 2 )\E^ 2 ) of <Sjv/2 ® $N/2> the latter symmetry implies 
that My = M\,. Diagonalizing My, one sees that 1^} can be 
written in the following form 

N/2+1 

|*>= £ Xi\w*)\wi) (24) 

with A; G R and \u>i) G Sjv/2 being eigenvalues and eigen- 
vectors of My, respectively. 

The fact that |^) G K(p Ts ) implies that for any pair 

\x), \y) G 5^/2, one has 

{x\y\p rs \^)=Y J Mx\y\p Ts H,u l ) 
1 

= y^ j Xi{wi,y\p\x,wi) 
1 

= A; {wi,y\p\ui, x) = 0, (25) 
1 

where the third equality follows from the fact that V s -g-p = p. 
As a result 

Tr[(W®\x){y\)p}=0 (26) 

holds for any pair of vectors \x), \y) G S^n, where W = 

On the other hand, with the aid of Eq. (24), one can rewrite 
Eq. (22) as 

(e\® N / 2 W\e}® N / 2 = 0. (27) 

Assume, in contrary, that the above (equivalently Eq. (22)) 
does not have any solution. Then, its left-hand side must have 
the same sign for any |e) G (D 2 , say positive. This means 
that the matrix W is an A^/2-qubit entanglement witness sup- 
ported on the symmetric subspace <Sjv/2. We already know 
that there are no PPT entangled symmetric states of two and 
three qubits, and hence for N = 4 and N = 6, this witness 
must be decomposable. Precisely 

W = P + Q Ta (28) 

for N = 4 with P, Q being positive matrices acting on (D 2 ® 
(D 2 , while 



( e |®*/2(e*|®w/2|^ =0 



(22) 



W = P + Q Ta + R Tab 



(29) 
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for N = 6 with P,Q,R>0 acting on (D 2 <g) (D 2 <g) (D 2 . Putting 
Eqs. (28) and (29) to Eq. (26), one in particular arrives at the 
following conditions 

Tr [(P ® |a?X^D p] = Tr [(Q ® \x)(x\) P Ta ] = 0. (30) 

for N = 4 and 

Tr[(P ® = Tr[(Q <g> \x){x\) P Ta ] 

= Tr[(R® \x){x\) P Tab ] = (31) 

for N — 6 and for any |a;) G <Sjv/2- These conditions imply 
that either p or p TA , or p TAB is not of full rank, contradicting 
the assumption. Therefore, Eq. (27) and thus Eq. (22) must 
have a solution. In this way we have proven the following 
theorem. 

Theorem 5. Four-qubit symmetric states of ranks (5,8,8) 
and six-qubit symmetric states of ranks (7, 12, 15, 15) are not 
edge. 

We already know that for there exist PPT entangled sym- 
metric states consisting of more than three qubits [8, 17, 18] 
and thus indecomposable entanglement witnesses detecting 
them. Consequently, the above method does not apply in gen- 
eral for even N > 8. Nevertheless, it provides a necessary 
condition for being edge, i.e., if a symmetric state of all ranks 
maximal except for r{p TN /' 2 ) which is one less than maximal 
is edge, the witness W = J^i l w z)( w i I constructed from (24) 
is indecomposable. 

Interestingly, for N = 4 we can prove analogous theorem 
in the case when the rank of r(p TAB ) is two less than maximal. 

Theorem 6. Generic four-qubit symmetric states of ranks 
(5, 8, 7) are not edge. 

Proof. By assumption, p and p TA are of full rank, while 
K(p TAB ) contains two linearly independent vectors (i = 
1,2). Consequently, to find a product vector e)® 4 G R(p) 
such that | e*)® 2 \ e)® 2 € R(p T - AB ), one has to solve two equa- 
tions 

(e*|® 2 (eH^) =0 (i = 1,2). (32) 

Let us now briefly characterize the vectors | With the 
aid of the identity p TAB = Vab,cd(p*) Tab Vab.cd, where 
Vab,cd is a unitary operator swapping AB and CD subsys- 
tems, one may show that they can be written as 

2 2 

|*i>=£A fc |e fe }|/*>, l*2)=E A ^>l e fc>- ( 33 ) 
fe=i fe=i 

Indeed, let us first notice that we can assume that one of 
is of Schmidt rank two. The largest subspace of (D 3 <g> (D 3 con- 
taining only vectors of Schmidt-rank three is one-dimensional 
(see, e.g., Ref. [34]). On the other hand, if one of |\I>j) (i = 
1, 2) is of rank one, i.e., is product with respect to the partition 
AB\CD, {e,r\p TAB \e,f*) = {e* J*\p\e* , f*} = mean- 
ing that p\e* , /*) = 0. It is then clear that since |e), |/) G <S 2 , 
P^e*,/*) G K(p), implying that r(p) = 4, which contra- 
dicts the assumption. 



Assuming then that is of rank two, either Vab.cd\^\) 
is linearly independent of leading to Eq. (33), or 

Vab,cd\^*) = for some £ G <D. In the latter case, 

short algebra implies that (i = 1, 2) are not linearly inde- 
pendent contradicting the fact that they span two-dimensional 
kernel of p T - AB . 

As a result, there is a vector \e) = (l,a) G (D 2 such that 
\e*)® 2 \e)® 2 G R(p TAB ) iff there is a G C solving the equa- 
tion 

P(a*)Q{a) + P{a*)Q{a) = 0, (34) 

where P, P and Q, Q are polynomials generically of degree 
two. Such a exists if and only if there is z G (D fulfilling 

P{a*)=zP{a*) (35) 

and 

Q{a) = -zQ(a). (36) 

With the aid of the first equation, we can determine a* as a 
function of z. There are clearly at most two such solutions. 
Putting them to Eq. (36) and getting rid of the square root, we 
arrive at a single equation 

{z*) 2 Q 4 (z) + z*Q' 4 (z) + Q'i{z) = 0, (37) 

where Q4, Q' 4 , and Q" stand for polynomials which are gener- 
ically of fourth degree. It has at least one solution because the 
polynomial appearing on the left-hand side of Eq. (37) has 
unequal degrees in z and z*. This allows for application of 
lemma 3, completing the proof. □ 

One notices that this method cannot be directly applied in 
the case of larger even N. Already for N = 6, the left-hand 
side of Eq. (34) contains three terms and therefore the factor- 
ization (35) and (36) cannot be done. Let us notice, however, 
that the numerical search for extremal states of even N done 
below does not reveal examples of extremal PPT entangled 
symmetric states of these ranks, suggesting the lack of edge 
states in generic states of these ranks. 

More generally, it should be noticed that the analysis of 
edge states allows for further reduction of configurations of 
ranks relevant for characterization of PPT entanglement in 
symmetric states. This is because a PPT state that is not edge 
can be written as a mixture of a pure product vector and an- 
other PPT state of lower ranks (see also Sec. V). 



C. On the Schmidt number of symmetric states 

Let us finally comment on the Schmidt number of the sym- 
metric states. Clearly, a pure state \ip) G T-Ln can be written 
as a linear combination of fully product vectors from T-Ln- 
Following Ref. [35], the smallest number of terms in such 
decompositions of is called the Schmidt rank of \ip) and 
denoted r(\ip}). Then, analogously to Ref. [36], we can de- 
fine the Schmidt number of p to be min/wA\{maXj r(\tjji))}, 
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where the minimum is taken over all decompositions {^i)} 
of p, i.e., p = J2i l^iX^tl- 

Below we show that in small symmetric systems consist- 
ing of four or five qubits, any entangled state has the Schmidt 
number two or at most three, respectively. We also comment 
on the Schmidt number of larger systems. 

Before that we need some preparation. Let us introduce the 
following transformations: F n : ((D 2 )®" -)• (D 2 and G n : 
((D 2 )®0 +1 ) h- ((Tj2)®(2n+i) defined through 



F„(l,a) K * = (l,a") 



and 



G n [(l,a n+i )® (1,q) 



(1,«) C 



(2n+l) 



(38) 



(39) 



respectively, for any a £ (D and n = 1,2, Notice that 

both maps are of full rank. Then, by F n and G n we denote 
maps that are defined through the adjoint actions of F n and 
G n , i.e., X(-) = X(-)Xt (X = F n ,G n ). 

Let us comment briefly on the properties of F n and G n . 
First, consider an A-qubit symmetric vector An ap- 

plication of -PYjv/21 to chosen [A/2] qubits of say the 
first ones, brings it to an ([A/2J + l)-qubit vector \ijj') £ 
(D 2 (g> <Si jv/2J ■ A subsequent application of Gjv/2-1 to the first 
A/2 qubits of \ip') in case of even N and Grjv/2i-i to the 
whole in the case of odd N, returns 

Analogously, by applying Fnv/21 to the first [A/2] qubits 
of an A-qubit mixed symmetric state p, one brings it to an 
(LA/2J + l)-qubit state a whose the last [A/2] qubits are 
still supported on the symmetric subspace «Si jv /2J ■ Let us de- 
note the parties of this state by Si, ... , B\j\r/2\+i- With re- 
spect to the bipartition B\ \B 2 . . . B^ N / 2 ^ +1 , it can be seen as 
a bipartite state acting on (D 2 ® (CL Ar / 2 J+ 1 . Moreover, this 
"compressing" operation preserves the rank of any symmetric 
state, i.e., r(er) = r(p). This is because R(p) is spanned by 
the symmetric product vectors (1, a)® N for a £ <D, which are 
then mapped by F { N/2] to (1, a)® W2J ® (l )Q; W21), Since 

all powers of a from the zeroth one to a still appear in the 
projected vectors, the whole information about p is encoded 
in a. Precisely, by an application of Gjv/2-1 to the first A/2 
qubits of a for even N and Gpjv/2]-i to all qubits for odd N 
of a, one recovers p. With all this we are now ready to state 
and prove the following theorem. 

Theorem 7. Let p be an entangled N-qubit symmetric state. 
If by an application of -Fnv/21 to the first [A/2] qubits of 
p one gets an ([A/2J + l)-qubit state a (see above) that is 
separable with respect to the bipartition Bi\B 2 ■ ■ ■ By N / 2 \+i> 
then p can be written as 



K 



[N/2] 



(40) 



where Aj, a* £ (D and [ip] denotes a projector onto 

Proof. We can clearly assume that r(p) = A + 1. Let the 
([A/2J + l)-qubit state a, coming from the application of 



F\n/2\ to the first [A/2] qubits of p, be separable with re- 
spect to the Hilbert space (D 2 <8 C LAr/2 J +1 . It can then be 
written as 



K 



Pi > 0, 2^pi = i, 

(41) 

where \e,) £ C 2 and |/,) £ CW2J+1 ^ 5 LAr/2J . 

As already said, the mapping ^rjv/2] preserves the rank, 
and therefore r(cr) = r(p) = N + 1. On the other hand, 
a acts on (D 2 (g> CD L JV / 2 J + 1 and therefore its maximal rank is 
2( [A/2] +1), which means that for odd A, a is of full rank, 
while for even A it has a single vector \<f>) = 1 10} — |0, A/2) 
in its kernel. 

Let us now divide our considerations into the cases of odd 
and even A. In the first case we put e^) = (1, a[ N ^ 2 ^ ), and 
then notice that the isomorphism (CL JV / 2 J+ 1 = S[n/ 2 ] im- 
plies that each |/j) £ (DL Ar / 2 J+ 1 can be expanded in terms 
of [N/2\ + 1 product symmetric vectors from iSnv72| m the 
following way 



\fi) 



[N/2\+l 



® LJV/2J 



with 



2nj 



and £ (D. Consequently, 



G\N/2-\-l\e-i)\fi) 



LJV/2J+1 



(42) 



(43) 



(44) 



with 



which, when substituted to Eq. (41) leads directly to (40). 
In the case of even A, a has a single vector \<f>) = |01) 

I A/2,0) in K(a). Putting again e. t ) = (l,af " 
aj £ (D, it imposes a constraint on that is, (N/2\fi) = 

(0\fi)a^ 2 , meaning that \ f) £ C^ 2 . Therefore, again ex- 
ploiting the isomorphism <D N ^ 2 = <Sjv/2-i> one sees that all 
can be written as 



7V/2 

i/ i )=E4 ) ( 1 ' e% «o w - 



(45) 



By substituting this to Eq. (41) and applying Gjv/2-1 to cr, 
one gets Eq. (40), completing the proof. 

Let us finally comment on the choice of the product sym- 
metric vectors used to expand |/;)'s. In both cases of odd 
and even A they are chosen to be (1, e'^'a)®^/ 2 ] (j = 
0, . . . , [A/2] — 1), where ipj are chosen so that c lipj are 
[A/2]th roots of the unity. One checks by hand that such 
vectors span a [A/2] -dimensional linear space. □ 

Theorem 7 implies that any PPT entangled symmetric state 
consisting of four (five) qubits has Schmidt number two (at 
most three). To be more precise, we prove the following corol- 
laries. 
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Corollary 1. Any PPT entangled symmetric state p of four 
qubits can be written as 

p=fyAf{l, a r i + 4\^-*iT% (46) 

i=l 

while any PPT entangled symmetric state of five qubits as 

i=l 

+4 t \l,e i ^a i r% (47) 
with K < 6 and ,tt,eC. 

Proof. By applying Fs (F3) to the first two (three) qubits of 
p for N = 4 (N — 5) we get a state a acting on C 2 ® (D 3 . 
It is clearly PPT and due to Ref. [30] also separable. Conse- 
quently, theorem 7 implies that p can be written as in Eq. (40), 
which in particular cases of A*" = 4 and N — 5 leads to (46) 
and (47), respectively. The number of elements in both the 
decompositions (46) and (47) follows from the fact that any 
qubit-qutrit separable state can be written as a convex combi- 
nation of six product vectors [33]. □ 

It should be noticed that by using the approach developed in 
Ref. [33], one can obtain decompositions of any PPT entan- 
gled four-qubit symmetric state similar to (46), but in which 
vectors (1, -a,)® 4 are replaced by either (0, l) m or (1, 0)® 4 . 

Theorem 8. Let p be a PPT entangled symmetric four-qubit 
state. Then it can be written as 

K 

p = Y,iMh «0® 4 + B t (p,in (48) 

i=l 

where K < 6, (l,a,) £ (D 2 , and Ai,Bi are some complex 
coefficients, and by [ip] we denote a projector onto 

Proof. The proof exploits the method developed in Ref. [33]. 
First, one notices that any p can be written as a sum of rank- 
one matrices 

K 

p = ^|^)(^|, (49) 

?:=i 

where, in particular, "1^) can be (unnormalized) eigenvectors 
of p, and K < 6 (see corollary 1). On the other hand, p can 
always be expressed in terms of the symmetric unnormalized 
basis {|-Ev)}„=i spanning £4 as 

5 

p= J2 P^\E*)(Et\- (50) 

Both decompositions (49) and (50) are related via the so- 
called Gram system of p, i.e., a collection of A'-dimensional 
vectors \v,) = (1/ (E*\ \E*), . . . , <* A -|£ 4 » (p = 
1, ... ,5), giving p^ = {vn\v v ). Putting the latter to Eq. (50) 
with explicit forms of the vectors \v^), one recovers (49). 



Now, by projecting the last party onto |0) we get a three- 
qubit symmetric PPT state p, which, as already stated, is 
separable. Then, according to Ref. [33], there exists a di- 
agonal matrix M = diag[a*, . . . , a* K \ such that \v^) = 
M^ 1 \v\) (p = 1, . . . , 4). For convenience we can also put 
U5) = M 4 \vi) + \v) with \v) being some A'-dimensional 
complex vector. Then, putting \vi) = (A\ : ... 1 A* K ) and 
\v) = (B*, . . . , B* K ), one sees that 

= A i J2ar 1 \Ef l )+B i \El) 

= A i {l, CH )®*+B i \Ef), (51) 

where the second equation follows from the explicit form of 
the vectors \v^). Substituting vectors to Eq. (49), one 
gets (48), which completes the proof. □ 

In order to get the representation (48) with (0, 1) replaced 
by (1, 0), one has to project the last party of p onto |1) instead 
of |0). 



IV. EXTREMAL PPT ENTANGLED SYMMETRIC STATES 

We will here seek extremal elements in the convex set of N- 
qubit symmetric PPT states Dpp^. We see that the number 
of distinct configurations of ranks for which one finds such 
examples is small and does not increase with N. In particular, 
if N is even there is only a single such configuration. 

Let us start by adapting to the multipartite case an algorithm 
described in Ref. [3] (see also. Ref. [7]) allowing to look for 
extremal elements of -Dppx- 



A. An algorithm allowing to seek multipartite extremal PPT 
states 

Let us consider again a product Hilbert space Hn = C dl (S> 
...(g) (D dl and the set -Dppt of fully PPT states acting on Hat. 
We call an element p of Dppt extremal iff it does not allow 
for the decomposition 

p = ppi + (1 - p)p 2 (52) 

with pi G -Dppt such that p\ ^ p 2 and < p < 1. Should 
p be not extremal, Eq. (52) implies that R(pi) C R(p) and 
R{P2) C R(p) and also R{pJ k ) C R(p Tk ) for all k. Alterna- 
tively speaking, if p is not extremal one is able to find another 
PPT state <t 7^ p such that 

R{a) C R(p) (53) 

and 

R(o T >) QR{ P Tk ) (k = l,...,M), (54) 
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where by M we denote the smallest number of partial trans- 
positions necessary to define the set of PPT states. 

On the other hand, if there is a PPT density matrix a ^ p 
fulfilling the conditions (53) and (54), then one finds such a 
number x > that p(x) = (1 + x)p — xa is again a PPT 
state. Thus, p cannot be extremal because it admits the de- 
composition (52), i.e., p = [p(x) + xa]/ (1 + x). In this way, 
one arrives at a natural criterion for extremality saying that p 
is extremal iff there is no state a ^ p such that (53) and (54) 
hold. 

Interestingly, one can even relax the assumption of a being 
a PPT state to being a Hermitian matrix satisfying a set of 
equations following from Eqs. (53) and (54): 

V k h Tk V k = h T « (fc = 0,...,M), (55) 

where To is an identity map. This is because if such h exists, 
one considers a one-parameter class of matrices p(x) = (1 + 
xTrh)p — xh (Trh is added in order to ensure that Tr[p(a;)] = 
1 for any x). It is straightforward to show that there exist 
x% < and X2 > such that for any x £ [xi,^], p{%i) 
(i = 1,2) are PPT states. Consequently, p = [l/(|aii| + 
\x2p{xi) + Xip(x2)], meaning that it cannot be extremal. 
Let us also notice that the system (55) is equivalent to a 
single equation 

[V M o...ot 1 oV]{h) = h, (56) 

where the maps Vk are defined through Vk{-) = 
[Pk{-) Tk P k ] Tk (k = 0, . . . , M). Indeed, if h satisfies the sys- 
tem (55) then it also obeys (56). On the other hand, if the con- 
dition (56) is satisfied with some h, then it has to obey (55); if 
one of the conditions (55) does not hold, a simple comparison 
of norms of both sides of (56) shows that (55) neither can be 
satisfied. 

As a result, we have just reached an operational criterion 
for extremality of elements of -Dppt [3] (see also Ref. [7]). 

Theorem 9. A given PPT state p acting on Hn is extremal 
in -Dppx iff there does not exists a Hermitian solution to the 
system (55) which is linearly independent of p. 

This also leads to a necessary condition for extremal- 
ity that can be formulated in terms of the ranks r(p Tk ) 
(fc = 1,...,M), Namely, each equation in (55) imposes 
[dim Hn] 2 — [ r {p Tk )] 2 linear constraints on the matrix h. The 
maximal number of the constraints imposed by the system is 

thus Efcio([ dim ^v] 2 _ \ r {p Tk )] 2 )- ° n the other hand > a 
Hermitian matrix acting on Hn is specified by [dim Hn] 2 real 
parameters and therefore if 

Yy{p Tk )] 2 > M[dh-nH N ] 2 + 1, (57) 

fc 

the system (55) has a solution and p is not extremal. 

Importantly, the above considerations imply an algorithm 
allowing to seek extremal elements of -Dppt [3]. Given a 
PPT state p and a solution h to the system (55), one considers 
p(x) = (1 + xTrh)p — xh. It is fairly easy to see that there 
is x = x„ for which p\ = p(x*) G -Dppt but r{p 1 k ) = 



r(p Tk ) — 1 for some fc. In other words, we can choose x 
in such way that one of the ranks of the resulting state p\ is 
diminished by one. 

For the resulting state p\ we again look for solutions to (55). 
If the only solution is p\ itself (up to normalization), then 
it is already extremal. If not, one again considers p\(x) = 
(1 + xTrh)px — xh and finds such x = x* that p 2 = Pi(x*) is 
a PPT state with one of the ranks diminished by one. We pro- 
ceed in this way until we obtain an extremal state. If the latter 
has rank one it is separable, otherwise it has to be entangled. 
Clearly, since we deal with finite-dimensional Hilbert space, 
the final state of this algorithm is reached in a finite number 
of steps. 

In our implementation of this algorithm we obtain the Her- 
mitian matrix h by solving a slightly different system of equa- 
tions than (55) (or the single one (56)), i.e., 

h T *\*f)=0 (i = l,...,k(p T *)) (58) 

for all fc = 0, . . . , M with ) denoting vectors spanning 
the kernel of p Tk . Clearly, both systems (55) and (58) are 
equivalent. 

A general Hermitian matrix acting on Hn is fully charac- 
terized by (dim "Hat) 2 real parameters hi, which we consider 
elements of a vector \h) G (C( dim «N) 2 . Then, the set (58) can 
be easily reformulated as a single matrix equation lZ\h) = 0, 

where TZ is a [2 dim 'Hat Y^lLq k(p Tl )] x (dim "Hat) 2 matrix 
with real entries. The number of rows of TZ stems from the fact 
that each group of equations in (58) corresponding to a given 
fc gives 2k(p Tk ) dimHN linear and real conditions for hi, 
which when summed over all partial transpositions (includ- 
ing also T ) results in the aforementioned number of rows. 
Clearly, this gives more conditions for hi than those follow- 
ing from Eqs. (55), i.e., ^ fc=0 [(dimH) 2 — [r(p^)] 2 ] (some of 
them are clearly linearly dependent). Nevertheless, even if the 
number of equations is larger in comparison to Eq. (56), our 
approach does not require multiplying many matrices, which 
in the case of many parties makes the implementation a bit 
faster. 



B. Extremal PPT entangled symmetric states 

Let us now apply the above considerations to the symmetric 
states. First of all, one notices that the condition (57) has to 
be slightly modified because different partial transpositions of 
a symmetric state and the state itself act on Hilbert spaces of 
different dimensions (see Ref. [7] for the case of four-qubits). 
Specifically, an equation in (55) corresponding to the partial 
transposition with respect to S gives (|5| + l) 2 (iV — \S\ + 
I) 2 - Hp Ts )? (\S\ = 0, . . . , [N/2\ ) linear equations. Then, 
the condition (57) becomes 

LA72J [N/2\ 

E t r ^ Ts )] 2 ^ E (|S|+i) 2 (AHS|+i) 2 -(vv+i) 2 +i. 

|S|=0 |S|=0 

(59) 
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Taking into account the fact that among symmetric states only 
those of rank N + 1 can be entangled, allows us to rewrite the 
above inequality as 

[N/2\ [N/2] 

E Kp Ts )] 2 > E (|s|+i) 2 (AHS|+i) 2 -(iv+i) 2 +i. 

|S|=1 |S|=1 

(60) 

In table I we collect all the ranks, singled out with the aid 
of Ineq. (60), for which there are no extremal PPT states for 
exemplary cases of low number of qubits iV = 4, 5, 6. 



N 


Inequality (60) 


Ranks excluded with (60) 


4 


[r(^)] 2 + [r(^)] 2 >121 


(5, 7, 9), (5, 8, 8), (5, 8, 9) 


5 


Hp Tl )] 2 + Hp T ' 2 )] 2 > 209 


(6, 9, 12), (6, 10, 11), (6, 10,12) 


6 


Wp Ti )] 2 + Kp T2 )] 2 

+[r(p T3 )] 2 > 577 


(7, 10, 15, 16), (7, 11, 15, 16), 
(7, 12, 14, 16), (7, 12,15,15), 
(7, 12, 15, 16) 



TABLE I. Inequality (60) (second column) for exemplary cases of 
N = 4, 5, 6 together with the ranks (third column) excluded with 
its aid for which there are no extremal entangled symmetric states. 
Notice that the fact that there are no extremal states of maximal ranks 
can be inferred without restoring to Ineq. (60). 



We have applied the above algorithm to the PPT symmet- 
ric states with the number of qubits varying from N = 4 to 
= 23. As the initial state we took the projector onto the 
symmetric space Vn- Clearly, it has all the ranks maximal, 
which is important from the point of view of the algorithm; 
if the ranks of the initial state p are too low (cf. theorems 2 
and 3), the algorithm cannot produce an entangled PPT state 
out of p by lowering its ranks. Then, we have searched for 
random extremal states by choosing randomly at each stage of 
the protocol the matrices h resulting from solving (55). On the 
other hand, to find other examples of ranks than those obtained 
through a random search and not excluded by the analysis 
above, we designed the matrices h in such way that they lower 
specific ranks. The obtained ranks are collected in table II. In- 
terestingly, there are always at most three different configura- 
tions of ranks assumed by the found extremal PPT entangled 
symmetric states and it seems that the number of configura- 
tions does not increase with N. Moreover, in the case of odd 
N, there is only a single such configuration (all ranks are max- 
imal except for the last one which is two less than maximal). 
It is an interesting problem to confirm these findings analyti- 
cally. And if this is the case, the problem of characterization 
of PPT entanglement in symmetric states reduces significantly 
to the characterization of extremal states assuming few differ- 
ent configurations of ranks, in particular a single one for odd 
N. Notice that in the case of symmetric qubits, there cannot 
be PPT entangled states of lower ranks than those assumed by 
extremal one as this is the case in higher-dimensional Hilbert 
spaces (cf. Ref. [4]). This is because in order to construct such 
states one needs PPT extremal entangled states supported on 
lower-dimensional Hilbert spaces, and in our case such states 
are always separable. 

Let us study in detail extremal entangled states in the exem- 



plary case of N — A. From theorems 2 and 3 it follows that 
PPT states of ranks (5, r(p TA ), r(p TAB )) with r(p TA ) < 6 or 
t(p Tab ) < 6 are either all separable or genetically separa- 
ble. Then, theorem 6 states that generic PPT states of ranks 
(5, 8, 7) are not edge and thus not extremal. Finally, Ineq. (56) 
implies that PPT states of ranks (5, 7, 9), (5, 8, 8) and (5, 8, 9) 
cannot be extremal. As a results, the natural candidates for 
extremal states that can be obtained with the aid of the above 
algorithm have ranks (5, 7, 7) and (5, 7, 8). 

We have run the algorithm 30000 times and 19.2% of the 
generated examples were extremal entangled states of ranks 
(5, 7, 8). In the remaining 80.8% of cases we arrived at states 
of ranks (5, 7, 7), all being separable. Also, when lower- 
ing the ranks from the initial state of ranks (5, 8, 9), 99.4% 
of the times we have obtained an intermediate (5, 8, 8)-state, 
whereas intermediate states of ranks (5,7,9) have appeared 
the 0.6% of remaining cases. Concluding, it should be no- 
ticed that with the aid of the above algorithm we have gen- 
erated PPT entangled extremal states assuming only a single 
configuration of ranks. All states of ranks (5, 7, 7) occurred to 
be separable and there is an indication suggesting that generic 
four-qubit PPT symmetric states of these ranks are separable 
(see appendix A). 

Let us finally notice that to make the application of the algo- 
rithm to systems consisting of even twenty-three qubits pos- 
sible, one has to take an advantage of the underlying symme- 
try and try to avoid representing a symmetric p and its partial 
transpositions in the full Hilbert space H2 N = ((C 2 )® N (see 
also Ref. [12]). Indeed, since Sn — < C N+1 , one can repre- 
sent p or a general Hermitian matrix supported on Sjy, as a 
(N + 1) x (N + 1) matrix, which we further denote p 1G d. In 
order to move from one representation to the other one we use 
a (N + 1) x 2 N matrix B N : S N i-> (D jv+1 given by 

N 

B N = E H(EZ\, (61) 
m=0 

which gives p re d = B^pB^. It is straightforward to check 
that for any N, B%B N = V N and B N Bj, = 1jv+i- 

Then, accordingly, the partial transposition of p with re- 
spect to Tfc can be represented as a (k + 1)(N — k + 1) x 
(k + 1)(N -k + 1) matrix p% d acting on <C k+1 ® € N ~ k+1 . 
To get the latter from p rcc j without restoring to the represen- 
tation of p in the full Hilbert space %2,n^ one can utilize a 
(k + l)(N-k+l)x(N+l)matnxB k = (B k ®B N _ k )B%, 
i.e., 

p£d = \( B k ® B N _ k )p{Bl ® B T N _ k )] T * 

= [{B k ® B N _ k )B T N p, cA B N (Bl ® B T N _ k )] T * 

= (B k p Icd Blf k . (62) 

Short algebra shows that the elements of B k (k = 
1, [N/2\) are given by 

with i = 0, . . . , k, j = 0, . . . , N - k, and n = 0, . . . , N. 
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Consequently, to get effectively a partial transposition of 
p, one maps a (N + 1) x (N + 1) matrix p re( j with Bk, and 
subsequently performs a simple partial transposition on the re- 
sulting bipartite matrix [cf. Eq. (62)]. Accordingly, one also 



transforms the Hermitian matrices h appearing in the system 
(58). Notice that this approach allows us to reduce the algo- 
rithm complexity from exponential to polynomial in TV both 
in time and memory. Precisely, the estimated time complexity 
of our approach amounts to 0(N 6 ). 



N 


r{p) 


r{p Vl ) 


r(p F2 ) 


r(/3) 




r(/s) 


r(p r °) 


r(p r7 ) 




r(/ 9 ) 




r(p Fl1 ) 


4 


5 


1 (-1) 


8 (-1) 




















5 


6 


10 


10 (-2) 




















6 


7 


12 


14 (-1) 
14 (-1) 


14 (-2) 
13 (-3) 


















7 


8 


14 


18 


18 (-2) 


















8 


9 


16 


21 


23 (-1) 
23 (-1) 


23 (-2) 
22 (-3) 
















9 


10 


18 


24 


28 


28 (-2) 
















10 


11 


20 


27 


32 


34 (-1) 

34 (-1) 

35 (+0) 


33 (-3) 

34 (-2) 
32 (-4) 














11 


12 


22 


30 


36 


40 


40 (-2) 














12 


13 


24 


33 


40 


45 


47 (-1) 

47 (-1) 

48 (+0) 


47 (-2) 
46 (-3) 
45 (-4) 












13 


14 


26 


36 


44 


50 


54 


54 (-2) 












14 


15 


28 


39 


48 


55 


60 


62 (-1) 

62 (-1) 

63 (+0) 


62 (-2) 
61 (-3) 
60 (-4) 










15 


16 


30 


42 


52 


60 


66 


70 


70 (-2) 










16 


17 


32 


45 


56 


65 


72 


77 


79 (-1) 

79 (-1) 

80 (+0) 


79 (-2) 
78 (-3) 
77 (-4) 








17 


18 


34 


48 


60 


70 


78 


84 


88 


88 (-2) 








18 


19 


36 


51 


64 


75 


84 


91 


96 


98 (-1) 

98 (-1) 

99 (+0) 


98 (-2) 
97 (-3) 
96 (-4) 






19 


20 


38 


54 


68 


80 


90 


98 


104 


108 


108 (-2) 






20 


21 


40 


57 


72 


85 


96 


105 


112 


117 


119 (-1) 

119 (-1) 

120 (+0) 


119 (-2) 
118 (-3) 
117 (-4) 




21 


22 


42 


60 


76 


90 


102 


112 


120 


126 


130 


130 (-2) 




22 


23 


44 


63 


80 


95 


108 


119 


128 


135 


140 


142 (-1) 

142 (-1) 

143 (+0) 


142 (-2) 
141 (-3) 
140 (-4) 


23 


24 


46 


66 


84 


100 


114 


126 


136 


144 


150 


154 


154 (-2) 



TABLE II. The ranks of extremal states found by using the algorithm described in Sec. IV A. The first column contains the number of qubits, 
while the next six columns the ranks of p and its partial transpositions r(p Ti ) (i = 1, . . . , 11). Notice that there are no PPT entangled states 
with less qubits than four [11] (cf. Sec. Ill A). The negative numbers in parentheses denote the difference between the given rank and its 
maximal value (the lack of parentheses means that the rank is maximal). For all N there are at most three possible configurations of ranks with 
extremal states, and, interestingly, in the case of odd N there is always only one such configuration. 



V. SPECIAL CASES N=4. It follows from theorem 2 that the four-qubit sym- 

Here we summarize the obtained results for particular sys- 
tems consisting of four, five, and six qubits. 
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metric states are separable if either r(p) < 4 or r(p TA ) < 4, 
or r(p T - AB ) < 3. Then, theorem 3 implies that if either 
r(p TA ) < 6 or r{p TAB ) < 6, generic symmetric p is sep- 
arable. This leaves only six configuration of ranks (out of 
all possible 72 assuming that r(p) = 5) among which one 
may seek PPT entangled symmetric states: (5, 7, 7), (5, 7, 8), 
(5, 7, 9), (5, 8, 7), (5, 8, 8), and (5, 8, 9). 

Passing to edgeness, it is known that states of ranks (5, 8, 9) 
cannot be edge. Then, it follows from theorem 5 that all states 
of ranks (5, 8, 8) are not edge, while from theorems 4 and 6 
that generic states of ranks (5, 7, 9) and (5, 8, 7) are not edge. 
These theorems also show that a typical PPT entangled state 
assuming one of the above four configurations of ranks can 
always be brought, by subtracting properly chosen symmetric 
fully product vector, to a PPT entangled state of ranks either 
(5,7,7) or (5,7,8). Interestingly, with the half-analytical- 
half-numerical method presented in Ref. [8] as well as the 
numerical algorithm described in Sec. IV B, we found solely 
examples of PPT entangled states of ranks (5,7,8) (due to 
Ineq. (60) PPT states of ranks (5, 7, 9), (5, 8, 8), and clearly 
(5, 8, 9) cannot be extremal). All the states of ranks (5, 7, 7) 
found with the above algorithm were separable. This together 
with the analytical considerations enclosed in appendix A sug- 
gest that that symmetric four-qubit states of these ranks are 
generically separable. Provided this is the case, the analy- 
sis of PPT entangled symmetric states of four-qubits could be 
reduced significantly to the characterization of states with a 
single configuration of ranks (5, 7, 8). 

N=5. In this case theorem 2 implies that five-qubit PPT 
symmetric states are separable if either r(p) < 5 or r(p TA ) < 
5, or r(p TAB ) < 4 are separable. Then, theorem 3 says that 
if either r(p T - A ) < 8 or r(p TAB ) < 9, they are generically 
separable. Similarly to the case of N = 4, this leaves six 
(out of 120 possible under the assumption that r(p) = 6) 
ranks for which typical PPT symmetric states need not be sep- 
arable: (6, 9, 10), (6, 9, 11), (6, 9, 12), (6, 10, 10), (6, 10, 11), 
and (6,10,12). 

With the aid of theorem 4, one sees that five-qubit PPT sym- 
metric states of ranks (6, 9, 12), (6, 10, 11), and (6, 10, 12) are 
generically not edge (notice that due to Ineq. 60 for the same 
ranks PPT states cannot be extremal). Hence, analysis of PPT 
entanglement in in this case reduces to three configurations of 
ranks (6, 9, 10), (6, 10, 10), and (6, 9, 11). Interestingly, only 
in the second case we found examples of extremal states with 
the above numerical algorithm (see table I). 

N=6. Let us finally consider the case of six qubits. The- 
orem 2 states that such PPT states are separable provided 
that either r(p) < 6 or r(p TA ) < 6, or r(p TAB ) < 5, or 
r(p TABC ) < 4. Moreover, theorem 3 implies that they are 
generically separable if either r{p TA ) < 10 or r (p TAB ) < 12, 
or r(p TABC ) < 12. The number of the remaining configu- 
rations of ranks among which one may seek PPT entangled 
states is then 24 (out of all possible 2880 when assumed that 
r(p) = 7), which is considerably larger than the correspond- 
ing numbers for N = 4, 5. 

Then, from theorem 5 it follows that six-qubit PPT states 
of ranks (7, 12, 15, 15) are not edge and theorem 4 says that 
generic PPT states of ranks (7, 12, 14, 16) and (7, 11, 15, 16) 



are also not edge. This, together with the fact that states of 
maximal ranks, i.e., (7, 12, 15, 16) are not edge, allows to re- 
duce the problem of characterization of six-qubit PPT states to 
still quite large number of 20 configurations. There are, never- 
theless, only two sets of ranks for which, using the algorithm 
from Sec. IV A, we found extremal PPT entangled states. 



VI. CONCLUSION 

Let us briefly summarize the obtained results. Our aim was 
to characterize PPT entanglement in symmetric states. We 
have made a significant step towards reaching this goal, yet the 
complete characterization for the general case remains open. 

First, we have derived simple separability criteria for PPT 
symmetric states in terms of their ranks, complementing the 
criterion stated in Ref. [11]. Interestingly, these criteria im- 
ply that for most of the possible configurations of ranks, PPT 
symmetric states are generically separable, and PPT entan- 
glement may appear only in a small fraction of cases, van- 
ishing for large number of parties. Putting r(p) = N + 1, 
there are precisely ([N/2\ +1)1 such configurations out of 
Af!(iV/2 + l) 2 (L w / 2 J" Ar / 2 )+ 1 all possible ones. For the exem- 
plary cases of four and five qubits this gives six different sets 
(out of, respectively, 72 and 120 all possible ones) of ranks for 
which typical PPT states need not be separable. 

Second, we have singled out some of the configurations 
of ranks for which PPT symmetric states are generically not 
edge, allowing for further reduction of relevant configurations 
of ranks. This is because if a PPT entangled state is not edge it 
can be decomposed as a convex combination of a pure product 
vector and a PPT symmetric state of lower ranks. From this 
point of view the relevant configurations of ranks are those 
than cannot be further reduced by subtracting a product vector 
from the state. Again, in the particular case of small systems 
consisting of four and five qubits, PPT states of higher ranks 
are generically not edge, lowering the number of the configu- 
rations to treatable two and three, respectively, for N = 4 and 
N = 5. 

Finally, with the aid of the algorithm proposed in Ref. [3], 
we have searched for extremal PPT symmetric states. We have 
investigated systems consisting of 4 up to 23 qubits and en- 
countered a clear pattern behind the configurations of ranks 
for which we have found examples of extremal states. In 
particular, for even N, except for the cases of N = 6 and 
A = 8, there are always three configurations (following the 
same pattern) of ranks. Interestingly, for odd N there is only 
a single such configuration, i.e., all ranks of the state and its 
partial transpositions are maximal except for the last one (the 
partial transposition with respect to half of the qubits), which 
amounts to two less than maximal. This is somehow in con- 
trary to the intuition that the number of different sets of ranks 
for which one finds extremal states should grow. On the other 
hand, it indicates that the problem of characterization of PPT 
entanglement in symmetric states could further be reduced to 
just few different types of states. 
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Appendix A: The case of (5, 7, 7) 

We consider here four-qubit PPT symmetric states of ranks 
(5, 7, 7) and provide a possible way of proving that generic 
states of these ranks are separable. For this purpose we use 
the approach developed in Ref. [33]. First, recall that any 
state p acting on TL can be written as a convex combination of 
rank-one projectors 



fc=l 



(Al) 



where the unnormalized vectors \ipk) are in general 
nonorthognal (a particular example of such decomposition is 
the eigendecomposition of p). 

Denoting by \ek) orthonormal vectors spanning H, we see 
that any element of p in this basis can be written as 

i 

<frWi) = ^(eilM**!^) = H«j> (A2) 
fc=i 

for any i, j = 1, . . . , dim ft, where the Z-dimensional vectors 
\vi) are defined as 




(i = 1,..., dim ft). (A3) 



Consequently, p is the so-called Gram matrix, i.e., the matrix 
of scalar products of a set of vectors called further the 

Gram system of p. A different decomposition in Eq. (Al) 
leads to a different Gram system and all Gram systems of p 
are related via unitary matrices (if extended to a properly large 
Hilbert space). 

Let us now consider a four-qubit PPT symmetric state p 
of ranks (5, 7, 7) and find its Gram system together with the 
Gram systems of p TA , p Ts , and p T - AB , starting from p TA . 
Since, by assumption, the latter is positive and r(p A ) = 7, 
it admits a decomposition as in Eq. (Al) (for instance the 
eigendecomposition) with seven rank-one components, i.e., 



/c=i 



(A4) 



Denoting by B the standard product basis {\i,j, k,l)} of 
((D 2 )® 4 and utilizing the fact that every |<J/ fc ) in (Al) belongs 
to (D 2 (&1S3, one finds that the Gram system of p TA with respect 
to B consists of eight seven-dimensional vectors \a), . . . , \d) 



and |tt), . . 


, \d) whose elements are giver 


by 




= (* fe |0000), 






b k 


= (tffc|0001) = 


(* fc |0010) = 


(**|0100), 


Ck 


= (* fe |0011) = 


(* fc |0101) = 


(**|0110), 


dk 


= (*fc|0111), 






CLk 


= (* fe |1000), 






b k 


= (¥ fc |1001) = 


(* fc |1010) = 


(**|1100), 


Ck 


= (* fc |1011) = 


(**|1101> = 


(*fc|H10), 


dk 


= (*fclllll>. 







(A5) 

Let us then introduce the following 7x4 matrices 

A = (\a),\b),\b),\c)), B = (\b),\c),\c),\d)), B = 

(|&), |c), |c), \d)), with columns 
\d). Then p T - A can be written as 



(\a),\b),\b),\c)), and C 
given by the vectors \a), . 



( A A^B A^B A*C \ 

B^A BjB BjB 

&A &B B^B fitc 

V &A &B &B &C ) 

(A 

fit ct) 



B 
B 



(A6) 



where X^Y (X, Y = A,B, C) denotes a 4 x 4 matrix con- 
sisting of scalar products of vectors defining X and Y. We 
will then symbolically denote p A = (A B B C). 

In the same way we can represent p. Since it is symmet- 
ric it admits the form p = (A' B' B' C) with A', B', and C 
constructed in the same way as A, B, etc. from the Gram 
system of p. Now, since both three-qubit matrices (0|p|0) 
and (0\p TA |0) arising by projecting the first qubit of p and 
p TA onto |0), are equal, they have the same Gram systems 
and therefore there is a unitary U such that A' = UA and 
B' = UB. Then, since by multiplying by a unitary opera- 
tor all the vectors of the Gramm system of p one gets another 
Gram system, we can always set A' = A and B' = B. 

Analogously, one sees that the matrices = 
(l\p TA |1) (arising by projecting the first qubit onto |1)), and 
therefore there is a unitary U such that B = UB and C = 
UC. In conclusion, we see that p and p Tji can be represented 
as 



{ABB C), 



P 



(A B UB UC). (Al) 



By taking partial transposition of p with respect to A and com- 
paring it with the above representation of p T A , one gets some 
conditions for U : 



where the vectors l^k) are subnormalized 



AWB = B^A, 
B^UB = CU, 



AWC = B^B, 
BWC = C^B. 



(A8) 
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The same reasoning allows us to represent p B and p AB 
in the following way 

p Tb ={A ub B UC), 
p Tab = (A UB VB VUC) (A9) 

with V being a unitary matrix and UB = VB. Again, com- 
parison of p TAB given by Eq. (A9) to the partial transposition 
of p Ts with respect to A, gives further conditions 

AWB = B^A, A^VUC = BlUB, 

b^uWb = c^Wa, b^Wvuc = &b. 

Now, having Gram systems of p and its partial transposi- 
tions, we introduce the following matrix 

7 

fe=l 

with | ^k) and | denoting decompositions of p Ts and p TAB 
[as in (A4)] corresponding to their Gram systems introduced 
above, and [)?/>)] standing for a projection onto \ip). In terms 
of the Gram systems (A9), Q assumes the form 

Q = (A UB B UC A UB VB VUC). (A12) 



By careful counting of the dimensions, one notices that Q is 
an unnormalized density matrix acting on (D ® (D 3 with re- 
spect to the bipartition aAB\CD, where a denotes the auxil- 
iary subsystem [cf. (All)], while ABCD stand for the sub- 
systems of p. Also, by definition Q is of rank seven, and 
hence according to Ref. [31] Q is separable with respect to 
this bipartition provided that it is supported on (D 7 ® (D 3 and 
QT a AB > o Although we cannot prove the former condition, 
it is clear that generic Q is supported on (D 7 Cg> (D 3 . Then, after 
direct algebra and with the aid of conditions (A8) and (A10), 
one sees that the latter condition, i.e., Q TaAB > holds if the 
following two equations 

B^VUC = C ] UB, C^U^VUC = C ] UC (A13) 

are obeyed. Still, exploiting the explicit forms of B and C, the 
above conditions can be further simplified, leading to a set of 
equations for scalar product of vectors composing the Gram 
system of p. Some of them, by virtue of the Eqs. (A8) and 
(A10) are immediately satisfied. Provided that the remain- 
ing equations also hold, one has that indeed Q TaAB > and 
generic Q is separable. 

Then, it is clear that by projecting the auxiliary qubit a 
of Q onto |0), one recovers p. So, if Q is separable across 
aAB\CD, then p is separable across AB\CD, implying that, 
due to the underlying symmetry, it is fully separable. 
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